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Abstract 

We directly measure by femtosecond time-resolved x-ray diffraction the symmetry coherent 
phonon excited in bismuth by strong optical excitation. The magnitude of the Eg mode observed 
is 0.2 pm peak-to-peak, compared against the 2.7 pm initial displacement of the fully-symmetric 
Aig mode. The much smaller motion of the Eg mode is a consequence of the short lifetime of the 
electronic states that drive the atomic motion. The experimentally measured magnitude of the 
Eg motion allows us to rule out a previously suggested scenario for explaining the dynamics in 
bismuth that relies on strong coupling between these modes. 
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I. INTRODUCTION 



Strong optical excitation of solids on time scales much smaller than typical vibrational 
periods often leads to coherent structural dynamics that is distinct from any dynamics 
exhibited in equilibrium. The best known example is that of elemental bismuth, where the 
electronic state distribution created by short pulse optical excitation leads to large-amplitude 
coherent vibrations that cannot be observed under conditions of thermal equilibriunPSI. 
Similar effects have been observed in related materials, such as telluriumP^. 

Coherent vibrations can in principle mediate phase transitions that involve a change in 
unit cell structure. In these cases the period of the vibrational mode sets a lower limit on 
the speed of the transition. The simplest example of this is the Peierls charge-density-wave 
transition in an idealized one-dimensional material with one atom per unit cell and one 
valence electron per aton#. At low temperatures, this structure is unstable with respect to 
the formation of a superlattice distortion with a periodicity of twice the spacing between the 
atoms. Sudden electronic excitation of the system in this low-temperature phase results in 
a coherent oscillation of the atoms as the interatomic potential energy surface suddenly re- 
laxes back toward the unmodulated, high-temperature structure. This is qualitatively what 
happens in the above-mentioned real systems like bismuth and tellurium. At sufficiently 
high levels of electronic excitation, models and some experiments have indicated that the 
interatomic potential can relax fully to the normal, unmodulated phase resulting in a fast 
phase transitiori^^. 

It has, however, long been recognized that short-pulse lasers can also induce large- 
amplitude coherent optical mode oscillations that break symmetry operations present in 
the initial stat^^. These kinds of motion are critical to understand in the context of in- 
ducing phase transitions that reduce the symmetry of a crystal. Bismuth is well-known 
as a model system for ultrafast structural dynamics, and as such it provides an excellent 
opportunity to study in detail the dynamics of non-fully-symmetric coherent phonon modes. 

The unit cell structure of bismuth under equilibrium conditions is shown in Fig. [Tj The 
primitive unit cell is rhombohedral with lattice constant = 4.7574 A and a = 57.09° at 
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room temperatur^I^. In cartesian coordinates, we can represent the primitive cell vectors as 

( 1 V3 l\ 
ai = I -2^h, — ^dh, gC I (Ij 

a2 = I ^ah, — —ah, -c I (2) 

as = ^0, ^ah, ^cj (3) 

where ah = 4.54675 A and c = 11.90291 A are the unit cell constants of the corresponding 
non-primitive hexagonal cell. The basis consists of two atoms, positioned along the diagonal 
of the unit cell at (0, 0, ±Cc), with C = 0.2334. 

In bismuth there exist three independent optical phonon modes at zero wave vector: 
one Aig mode and a pair of degenerate Eg modes. The Ai^ mode preserves the symmetry 
of the crystal and can be described as an oscillatory movement of the two basis atoms in 
opposite directions along the body diagonal of the unit cell with an amplitude Uz- The 
eigenvectors for the modes are sketched in Fig. [ij these modes correspond to motion of 
the basis atoms in equal and opposite directions perpendicular to the C3 axis. Displacement 
along any coordinate breaks the C3 symmetry of the cell. A displacement Ux along the 
direction of the x-axis will also break all C2 and mirror plane symmetries. For the orthogonal 
Uy displacement, exactly one C2 axis and one mirror plane symmetry are preserved. 

Impulsive stimulated Raman scattering from an optical pulse with a duration much 
smaller than the period of these modes can create large amplitude coherent motions of 
zero wave vector phonon mode^. Stevens et al.'^ proposed a treatment of this process for 
absorbing media that involves for each phonon mode two frequency-dependent second-rank 
tensors. One of these tensors that we will call vr^^ (j = x, y, z) is of interest here since it 
gives the force on an atom along the Uj directions: 

Fjit) = ^-^J2 [ I e-'^'Ekiuj)TTi{iu,uj~ n)E*{uj - n)dudn + c.c. (4) 

kl 

where Ek{uj) is the Fourier transform of the optical electric field, and Vc is the volume of 
a unit celP^. Each phonon mode is associated with a different tensor vr^^. The symmetry 
properties of the phonon modes place restrictions on the form of vr^^. For the Aig mode 
the symmetry of the crystal is preserved for any value of the displacement u^, so the tensor 



takes the diagonal form 
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For the E„ modes, the C3 symmetry imphes there are two tensors 
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If the optical photon energy is far from electronic resonances the values of d and / are in 
principle measurable from spontaneous Raman scattering. Near resonances the situation 
becomes more complex since the dynamics of the polarization density rely on another non- 
equivalent tensor with the same symmetry form but different frequency dependenc^^. 

Coherent excitation of the Aig and Eg modes in bismuth has been previously observed 
with pump-probe optical reflectivity experiment j^ * ^^ * ^^ . One of the more intriguing outcomes 
of this work is the suggestion of coupling between the two mode^^'^. The strength of 
this coupling is highly sensitive to the amplitude of the coherent Eg mode. Consequently, 
evaluating possible coupling mechanisms requires some way to accurately experimentally 
measure the actual atomic displacements of both the Eg and Aig modes, quantities that 
are not available from transient reflectivity measurements. More generally, quantitative 
information on the magnitude of non-fully-symmetric modes permits investigation into the 
possibility of inducing structural phase transitions where the target has a symmetry lower 
than the initial structure. This requires a more direct type of measurement, now possible 
using femtosecond time-resolved diffraction technique j^ * ^ * ^^ l 



II. EXPERIMENT 

Figure [2] shows a conceptual sketch of the femtosecond x-ray diffraction measurements. 
An intense infrared pump pulse (800 nm, 115 fs, 1 kHz) generated by an amplified fem- 
tosecond laser system excites the sample. To probe the changes in structure, we use the 
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140 fs femtosecond duration pulses generated by the electron-beam slicing method at the 
Swiss Light Sourc^^. As in our previous investigations on bismutli ^^ * ^^ * ^^ *, a toroidally bent 
mirror first coUimates the beam vertically and focuses it horizontally to achieve a spot size 
of 250 at the sample. Further downstream, an elliptically bent mirror focuses the beam 
vertically to a size of 7 /im at the sample position. Before hitting the sample, the x-rays 
reflect horizontally from a M0/B4C multilayer that weakly monochromatizes the x-ray beam 
to 7.05 keV with a bandwidth of 1.2%. The x-rays enter the sample at a grazing angle of 0.5° 
in order to limit the x-ray field penetration depth to 30 nm and to better match the pump 
excitation depth. The surface of the sample was miscut from the (111) plane by = 57° 
toward the (2 1 1) plane. Rotating the sample about its surface normal allows access to both 
the (1 1 1) and (0 1) x-ray reflections. 

Previous femtosecond x-ray diffraction experiments on bismuth have used a pump beam 
with a large incidence angle in an effort to minimize the angle between the pump beam and 
the x-rays. The non-zero size of the x-ray beam projected onto the crystal surface causes 
a loss in effective time resolution from the difference in relative arrival times between the 
pump and x-ray beams at different points along the sample. For a 7 fim vertical focus size 
of the x-rays at a grazing angle of 0.5°, the time resolution from this geometrical factor runs 
from 40 fs at 80° pump incidence angle (10° grazing) to 2.6 ps at normal pump incidence. 
High incidence angles for the pump beam come, however, with a price: the transmission 
and reflectivity, especially near the Brewster angle, become highly polarization dependent. 
In the case of bismuth at 80° incidence, the reflectivity of 800 nm light changes from 33% 
for p-polarization to 94% for s-polarization. This makes it difficult to study the effects of 
changes in polarization alone on the dynamics, which is one of the key differences between 
the excitation of the Ai^ and modes. For this it is advantageous to make the pump beam 
arrive at normal incidence where the difference in transmission and reflectivity for different 
polarizations is minimized. 

To achieve this without compromising time resolution, we adopted a scheme to tilt the 
intensity front of the pump laser pulse by 45° to make the pump-probe delay time nearly 
constant for normal incidence across the entire pumped region of the crystal. Practically, 
this is done by imaging the first order reflection from a grating onto the sampld2Sl_ xhe tilt 
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angle /5 is set according to 

tan(/3) = , ^° = (8) 

dgM^Jl- {\Q/dg-sinQf 

where Aq is the center wavelength of the laser pulse, dg is the grating line spacing, B is the 
grating incidence angle, and M is the magnification of the imaging system. Both and M 
may be adjusted to tune (3. In our case we user a blazed gold-coated grating (Spectrogon 
PC0900) with dg = 900 lines/mm and 6 = 17°, and so to achieve (3 = 45° we adjusted 
the path delays to make M = 0.80 while maintaining the image point at the sample. The 
imaging was performed with a pair of plano-convex lenses (focal lengths /a = 1.5 m and 
/4 = Im) configured to compensate partially for spherical aberration^^. The alignment and 
tilt angle was verified prior to the x-ray experiment by cross-correlating an untilted pick-off 
of the 800 nm beam with the tilted beam using a nonlinear BBO crystal cut for phase- 
matched second harmonic generation at the optical position of the sample. The size of the 
laser spot on the grating was controlled by a pair of lenses /i and /2 that reduce the size of 
the beam to a diameter of 3 mm (FWHM). The incident fluence on the sample was set to 
6 mJ/cm^, corresponding to 1.6 mJ/cm^ absorbed fluence. 



III. RESULTS 

Because the (111) atomic planes lie perpendicular to the diagonal of the unit cell, x-ray 
diffraction from these planes is insensitive to atomic motion along the or Uy coordinates. 
The intensity of this reflection is, however, strongly sensitive to coherent motion of the Uz 
coordinat(™2ll3]_ xhe (0 1) planes lie at an angle of 71.6° from the (111) planes and are 
strongly sensitive to both the Ai^ mode and to the Eg mode with displacement along the 
Uy coordinate. Quantitatively, the intensity / of a diffraction peak with reciprocal lattice 
vector G in a kinematic approximation is proportional to where F = /je**^'""^ is the 
structure factor. Here the sum is over all basis atoms with index j, fj are the atomic form 
factors, and are the basis atom positions. For the speciflc case of bismuth and the (111) 
and (001) reflections we may write 



and 




cos2(27r( 



(10) 



cos2(27rC) 



where and /qqi are the measured intensities in equihbrium ( = 0). Here we 



lead to changes of less than 1% for the diffraction peaks we study here. 

Figure |3] shows the diffracted intensity from the (111) and (001) diffraction peaks as a 
function of pump-probe delay at an initial sample temperature of 170 K. For these mea- 
surements a half-wave plate controlled the polarization 6 of the pump with respect to the 
projection of the C3 axis with an uncertainty of ±2°. Positive values of 6 denote counter- 
clockwise rotation. In Figs.jS^a) and[3|d), 6 = —30°. In Figs.[3|b) and|3](e), this was changed 
to = 60°. There is little difference in the time-dependence of the (111) diffraction, but the 
dynamics from the (001) peak show noticeable changes. This becomes clear in Figs. |3](c) 
and [sj^f) which show the difference between the data sets when changing only the pump 
polarization. For the (111) diffraction the difference is zero within the errors, while for the 
(001) diffraction there is a clear oscillation starting at a pump-probe delay of zero. To fur- 
ther establish that these oscillations depend sensitively on the polarization. Fig. |4] shows for 
each diffraction peak the polarization dependence of the ratio of the pump-induced intensity 
change AI{ti)/ AI(t2), where ti = 80 fs and ^2 = 370 fs are pump-probe delays corresponding 
to the approximate times of the first two extrema of the oscillation. For the (111) peak this 
is independent of polarization within the experimental errors, whereas for the (001) peak 
there is a strong dependence with a period of 180°. 



IV. DISCUSSION 

A quantitative analysis of the polarization dependence requires us to evaluate the optical 
field inside the crystal, which is in general straightforward but algebraically complicated 
due to birefringence of the crystal. Fortunately, at 800 nm in bismuth the birefringence is 
fairly weak: en = £22 = —18.4 + 28^ and €33 = —13.5 + 28t^. We accordingly make the 
approximation of an isotropic dielectric tensor with e ~ ejj/S = —16.8 + 28i. If -E'o(^) is 



neglect treatment of changes in the Debye- Waller factor that, while present in bismuthl^^l^ 
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the amplitude of the incident electric field, inside the crystal we then have 

2Eo{t) 



^ ^ cos Om COS 6* — ^ sin 6* ^ 



1 + 



2 

sin am cos 9 



\ cos a-m cos 6^ + ^ sin 6^ 



(11) 



y sm am cos a j 

To determine the dependence of the phonon motion on the polarization angle 9, we 
can use Eq. |4] which gives the time-dependent force on an atom along one of the phonon 
eigenvectors. The resulting motion is then given by the classical equation of motion 



+^^>, + 27,^ = F,(t)/m (12) 

where m is the mass of an atom, Vtj is the mode frequency, 7j is the damping rate, and 
Uj is the position of the atom along the eigenvector. The initial conditions at t = — oo are 

du-i r\ 

= ^ = 0. 

We first consider the Aig mode, where the relevant coordinate is Uz- The behavior of 
this mode is usually described using the "displacive excitation" modeP^'^Sl^ ^he idea here is 
that electronic excitation causes a sudden displacement of the quasi-equilibrium value of the 
phonon coordinate away from zero, which in Eq. [12] is equivalent to a step-like behavior in 
Fz{t). If the lifetime of the electronic relaxation is much longer than that of a phonon period. 



the solution of Eq. 12 for is a cosine-like oscillation about the displaced quasi-equilibrium 
value of the coordinate. The dependence of the Uz motion on the pump polarization in the 
displacive excitation model is given simply by the dependence of the absorption of the pump 
light on the polarization. As noted above, in bismuth at 800 nm this dependence is relatively 
weak and leads to overall variation in the pump-induced displacement of approximately 10%. 

We now consider the Eg modes. In this case the modes are driven by electronic excitation 
to states that break some symmetry operations of the crystal, a subset of the states that form 
the conduction bands of bismuth. The lifetime of these states is then a critical parameter for 
characterizing the structural response of the system. If these states are long-lived, Stevens 
et al. argue that a similar argument would apply as for the Aig mode and yield a displacive 
form of the force on the Eg coordinates and a cosine-like displaced oscillation. This is, 
however, clearly not the case in the experiment: the Uy coordinate oscillates about zero, and 
has a phase closer to that of a damped sine function. Riffe and Sabbah'^ have proposed an 
extension of the model of Stevens et al. that can cover the case where the relevant electronic 
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state lifetime is very small, as might be expected from strong carrier-carrier scattering. The 
effective force on the Eg coordinate is then approximately proportional to the instantaneous 
intensity of the pump pulse inside the crystal. In the frequency-space formulation of the 
force given in Eq. |4| this is equivalent to approximating the elements d and / of '7Tjl}'(u,u — il) 
as independent of Q for values of Q comparable to or smaller than the inverse pump pulse 
duration. Applying Eq. [4] for the force along the Uy coordinate and solving for the motion 
of Uy gives 

/oo 
g{t - t')h{t')dt', (13) 
'OO 



where 



is the area-normalized intensity envelope of the laser pulse, 



git) = ^e-*^/-^ (14) 



,0 t<0 

h{t) = { , ^ , (15) 

4^6-^^* sin - -f^t t > 



C = A-4d, (17) 
D = VA^ + 52, (18) 
S = tan-i , (19) 



A = rf(3 + cos2a^) + 2/sin2a^, (20) 
B = AV3{—d cos am + f sin am) (21) 

and J-" is the incident laser fiuence. 

The amplitude of the displacement is the sum of a term that is polarization independent 
and a term that varies as cos(20 — 6). The polarization dependence of the time-resolved 
diffraction from (001) shown in Fig. |4[b) shows a cos(2^^ — 5) dependence in A/(ti)/A/(t2)- 



The solid line here shows a fit to the form of Eq. 16 The fitted value of 5 = — 1.04±0.14 gives 
us directly an estimate for the ratio f /d = — 0.07±0.14. Figure [s] shows the difference in the 
Uy dynamics between the polarizations 6 = —30° and 6 = 60°, inferred from the diffraction 



data using Eq. 10 The solid lines in Fig. p^are fits to Auy{t) = Uy{t, 6 = —30°) — Uy{t, 6 



60°) using the above expressions. Table 1 summarizes the resulting fit parameters. 



Our estimates of S and D allow us to estimate absolute values for the Raman tensor 
components d and / under the assumption that 6 is independent of temperature. These 
are listed in table 1. These values imply that the order of magnitude of 10~^ for changes 
in the dielectric constant due to the phonon motion, which is roughly consistent with the 
magnitude of reflectivity changes observed in optical reflectivity measurement^^*^^*^. The 
large decrease in signal when increasing temperature from 170 K to 300 K appears to be due 
to a large increase in the damping of the mode, suggesting that the mechanism of damping is 
related to interactions with incoherent phonon modes that are more populated at the higher 
temperature. 

We also note that the peak-to-peak amplitude of the oscillations (0.2 pm) is more than 
a factor of 10 smaller than the atomic motions associated with the Aig mode (2.7 ±0.1 pm) 
at the measured excitation fiuence. The difference may be understood by considering once 
again the lifetime of the electronic excitations that drive each mode. We can estimate the 
relative force on the excited phonon by taking the Fourier transform of Eq. |4]and evaluating 
its magnitude at the phonon frequency. Assuming states with a linear decay constant 7 
contribute to the phonon force is (for short driving pulses) 

For the Ai^ mode 7 ^ fi^, and so the amplitude of the phonon is approximately proportional 
to 1/^2. The short-lived electronic states that drive the E^ modes, however, have a much 
higher damping 7 ^ ^x,y The amplitude of the phonon is then proportional to I/7. If 
we assume the proportionality constants are of the same order of magnitude, from the 
relative size of the phonon amplitudes we can estimate I/7 ~ 4fs. This is consistent 
with our expectations based on the sinusoidal phase of the Eg oscillations. Note that the 
sinusoidal phase would also be expected for impulsive Raman scattering from coupling to 
virtual electronic transitions, which gives an additional purely real contribution to F{Q). 
The spectral dependence of the relative Raman cross section indicates, however, that the 
resonant contributions dominate in this range of excitation wavelengths^ . 

The measured magnitude of the coherent Eg modes allows us to test proposed models 
of phonon-phonon coupling in bismuth. In Ref. [3l Hase et al. report transient reflectivity 
measurements on bismuth that show possible evidence of coupling between the Ai^ and 
Eg modes at excitation levels approximately equivalent to those used in our diffraction 
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experiment' . In this work a low-frequency side-peak in the spectral density of the transient 
reflectivity appears with an amplitude of approximately 6% that of the main Ai^ peak. 
In Ref. 1201 Zijlstra et al. propose to explain this side-peak as a consequence of coupling 
between the Aig and Eg modes. They performed numerical simulations to support this 
idea. The simulated data reproduce a side-peak at a similar frequency as seen in the optical 
reflectivity data, but the magnitude of the peak depends strongly on the assumed coherent 
Eg amplitude. For a peak-to-peak E^ amplitude of O.lSAu^, where Au^ is the initial peak- 
to-peak amplitude of the Ai^ mode, the phonon-phonon coupling model predicts a spectral 
side-peak with a magnitude only 7x 10^^ that of the dominant Aig peak. Using our measured 
An^ = 2.7pm, we see that the peak-to-peak E^ amplitude is in reality COTSAu^, a factor 
of two smaller. We conclude that phonon-phonon coupling as presented in Ref. [20] is not 
sufficient to explain this feature of the optical reflectivity data. 

These results also confirm the interpretation of Ref. |2U where the oscillatory dynamics 
of diffraction from the (101) and (112) planes in bismuth were ascribed to excitation- 
induced changes in the Debye- Waller factor. Contributions from the E^ modes were excluded 
based on the time dependence of the observed dynamics, which were inconsistent with the 
expectations from coherent E^ oscillations. The present work allows us to also make an 
upper bound on the magnitude of coherent E^ contributions of 8 x 10^^ for the (1 1 0) peak 
and 3 x 10~^ for the (1 1 2) peak. These are much smaller than the 10~^ noise level on 
these measurements, implying that the observed dynamics are indeed related to changes 
in the Debye- Waller factor and not coherent E^ phonons. A small contribution of the 
approximately 1.4 THz Debye- Waller oscillations may also explain the low- frequency side- 
peak seen in the optical data of ref. |3] Alternatively, the side-peak could also be explained 
as a direct contribution to the reflectivity from the coherent E^ mode because of a small 
optical alignment error. 

As discussed in the introduction, coherent control over non-fully-symmetric optical 
phonons may be a viable route to the control of a variety of phase transitions. For this 
purpose much larger amplitudes of coherent excitation are desirable. One key to increasing 
the amplitude of the coherent motion under the resonant Raman scheme is to increase 
the lifetime of the electronic states that drive the transition. Greater amplitudes may be 
reahzed in semiconducting or insulating systems pumped just above bandgap, or in strongly 
correlated electron systems where the flnal electronic states are either long-lived or take sig- 
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nificantly longer times to equilibrate. Alternatively, direct dipole excitation using coherent 
mid-infrared or terahertz frequency pulses show promise in this area for cases where the 
mode is infrared active. 
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FIG. 1. Sketch of the rhombohedral unit cell of bismuth, showing the displacements corresponding 
to each zone-center optical phonon. The Ux displacement moves the atoms parallel to a C2 sym- 
metry axis, whereas the Uy displacement moves the atoms within a mirror plane containing the C3 
symmetry axis. 
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FIG. 2. Simplified sketch of the experimental setup. To excite the dynamics in the sample, a 
femtosecond pulse from a commercial regenerative amplifier is tilted by imaging the first order 
reflection from a grating onto the sample. The polarization of the pump is controlled by a A/2 
plate inserted just before the sample. To probe the dynamics, the sliced femtosecond duration 
x-rays reflect from an elliptically bent mirror vertically, and then horizontally from a multilayer 
mirror before encountering the sample. Diffracted x-rays are then detected using an avalanche 
photodiode. The reflection conditions for different diffraction peaks are realized by adjusting the 
sample rotation (j). 
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FIG. 3. Summary of time-resolved diffraction data from a bismuth single crystal at a temperature 
of 170 K. The left panels show data for the (111) reflection, whereas the right panels show data 
for the (100) reflection. The topmost panels show the data for a pump polarization = —30°; the 
middle panels show the data for 6 = 60°. The bottom panels show the difference in the diffracted 
intensity between the two polarizations. Whereas the (111) peak shows no significant difference 
between the polarizations, the (001) data shows evidence of a strongly damped oscillation. The 
arrows indicate the times that are used in Fig. |4] to study the polarization dependence in more 
detail. 
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FIG. 4. Polarization dependence of the ratio of the pump-induced changes between times ti = 80 fs 
and t2 = 370 fs after excitation, for (a) the (111) reflection and (b) the (001) reflection. While 
diffraction from (111) shows no dependence on the polarization, that from the (001) planes shows 
a strong dependence with a period of 180°. The solid curve shows a fit to a model discussed in the 
text. 
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FIG. 5. Difference in the Uy coordinate as a function of pump-probe delay time for excitation po- 
larizations at 6 = —30° and 6 = 60°. The different panels show data taken at sample temperatures 
of (a) 100 K, (b) 170 K and (c) 300 K. The solid lines are fits to a damped coherent Eg mode 
discussed in the text. 
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Temperature (K) Vy (THz) (ps ^) D (pm d (pm f (pm ^) 

100 1.84 ±0.06 2.9 ±0.5 10.5 ± 1.3 2.2 ± 0.4 -0.2 ± 0.3 

170 1.8 ±0.04 1.8 ±0.3 8.4 ± 1.2 1.7 ±0.3 -0.1 ±0.3 

300 1.9 ±0.7 12 ±10 9±7 1.8 ±1.1 -0.1 ± 1.0 

TABLE I. Summary of fit parameters for the data in Fig. [5} using the fit function discussed in 
the text. Here Vy = VlyjlTx. The estimated uncertainties reported for D, d and / do not take into 
account an additional 10% uncertainty in the excitation fluence, stemming from the uncertainty in 
the laser spot size that was the same for all measurements. 
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